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Введение. Гидромуфты (ГМ) нашли широкое применение в силовых приводах ленточных, цепных 
скребковых и пластинчатых конвейеров, элеваторов, осевых вентиляторов и дымососов, в приво- 
дах питательных насосов и газовых турбин, дробилок и мельниц различных типов, роторных экс- 
каваторов, дорожных катков, бетоносмесителей, барабанных сушилок и центрифуг, в составе 
гидромеханических коробок передач тепловозов, автомобилей, тракторов. Использование ГМ в 
составе силового привода даёт следующие преимущества: 

® улучшение статических и пусковых характеристик машины; 

® защита привода от динамических перегрузок при резкой остановке, например за- 
клинивании рабочего органа машины; 

® ограничение крутящего момента заданной величиной в режимах пуска и торможе- 
ния, что предохраняет двигатель, привод и рабочий орган машины от недопустимых перегрузок; 

® демпфирование и снижение амплитуды крутильных колебаний в приводах, возни- 
кающих при работе машины. 

Наиболее распространён электропривод с асинхронными двигателями (АД) [1]. В ряде 
случаев ГМ является здесь связующим звеном между силовым валом АД и приводимого механиз- 
ма. При этом ГМ обеспечивает плавность разгона механизма передачи энергии и облегчение ус- 
ловий прямого пуска АД на холостом ходу и с нагрузкой. 

Эффективным способом анализа динамических процессов, поиска путей совершенствова- 
ния конструкций является компьютерное моделирование. В связи с постоянным совершенствова- 
нием аппаратной части ЭВМ и программного обеспечения не перестаёт быть актуальной задача 
создания современных инструментов имитационного моделирования в электро- и гидроприводе, в 
том числе с использованием технологий визуального программирования. 

Математическая модель трёхфазного АД. В основу компьютерной модели АД с ГМ положена 
математическая модель трёхфазной электрической машины в заторможенной системе координат 


а, В, у [2], оси которой совмещены с осями А, В, Сфаз статора. Модель дополнена с целью учёта 
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потерь энергии в стали активными сопротивлениями г, контура намагничивания, включёнными в 
каждой фазе параллельно главным индуктивностям [, [3]. Уравнения, описывающие модель элек- 


тромагнитных процессов в АД, представлены ниже. 
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Здесь приняты обычные для математической модели обобщённой электрической машины 
допущения [2]: каждая из трёх обмоток статора создаёт в гладком зазоре синусоидально распре- 
делённую МДС. Коэффициент насыщения магнитной цепи принят неизменным. В выражениях (1) 
и далее приняты следующие обозначения: и — напряжение, / — ток, Ё — время, г — активное со- 
противление, Ф — потокосцепление, и — частота вращения. Нижние индексы а, В, у обозначают 
принадлежность к соответствующей фазе. Нижний индекс $ обозначает принадлежность к стато- 
ру, индекс г — принадлежность к ротору, а индекс и — принадлежность к ветви намагничивания 
Т-образной схемы замещения АД. /, асме — активная составляющая тока намагничивания. Детали- 
зация составляющих системы (1) приведена в выражениях (2) — (5). 
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Потокосцепления [1, 2]: 
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В выражениях (5) / — реактивная (индуктивная) составляющая тока намагничива- 


н_ геасвие 
ния фазы АД. Для АД типа АЖВ250М2ЗУХЛ2 мощностью 110 кВт активная составляющая тока на- 


магничивания / 


н_ геасвуе 


составляет приблизительно 5 % от тока намагничивания /.. Это справедли- 


во как в ходе переходного процесса пуска, так и при установившемся режиме. 


Уравнение электромагнитного момента АД имеет вид: 


и рее 
Ми =р` | 54 Г(4, йа (в +1): (6) 
Кинематическая схема механической части привода принята следующая: вал АД — повы- 
шающий редуктор — вал насоса — ГМ — вал турбины, к которому присоединена механическая 


нагрузка. Уравнение движения для вала АД имеет вид: 


4, р а ГМ Р=нехо 0, ‚ К 
СЁ Е. И] Мы 7 а © сопр. н®, |. (7) 
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Здесь Л — момент инерции вращающихся частей ротора, редуктора и насоса с рабочей 
жидкостью; р — число пар полюсов АД; передаточное отношение ], редуктора — это отношение 


частоты вращения ротора АД к частоте вращения насоса ГМ 


Тв = ©, О р (8) 
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М ососа гм |» — вращающий момент на валу привода насоса ГМ, приведённый к валу ротора АД; 


Риехо 


17 
[0 . 
|=) — момент механических потерь АД при способе вентиляции электрической машины 
о 


о 


«самовентиляция» [4]; р — механические потери АД в режиме идеального холостого хода; 


мех 0 


©&„, — частота вращения АД в режиме идеального холостого хода; К (©, — момент сил сопро- 


сопр. Нк 
тивления вращению валов насоса ГМ и редуктора, приведённый к валу ротора АД. 

Зачастую необходимо рассчитать энергетические показатели электрической машины в пе- 
реходном процессе. Для АД при отсутствии электрического питания обмоток ротора мгновенное 


значение активной потребляемой мощности Р вычисляется с помощью равенства 
Я =Иа + И И, (9) 


Мгновенное значение реактивной потребляемой мощности О, равно: 


о - [и (ыы) чи (1-1) /В. (10) 
Мгновенное значение полезной мощности Р, на валу АД: 


Р, м" И Й ; (11) 

Форма уравнений математической модели АД, представленная в трёхфазных координатах 
а, В, у, удобна при анализе и сравнении расчётных и реальных значений токов и напряжений фаз 
АД. Эта форма без дополнительных преобразований и уравнений пригодна для анализа особен- 
ностей несимметричных по характеристикам питания или параметрам режимов работы АД. 
Компьютерная модель АД. При создании компьютерной модели АД средствами САПР, предна- 
значенной в основном для проектирования и имитационного схемотехнического моделирования 
электронных и электротехнических устройств, была поставлена задача обеспечения встраиваемо- 
сти модели АД в единую комплексную модель электрогидродинамической передачи энергии. По- 
лучена компьютерная модель, пригодная для описания АД как с короткозамкнутым, так и с фаз- 
ным ротором. Обмотки фаз могут соединяться по схемам «треугольник», «звезда», иметь любую 
другую схему соединений между собой или с источником питания. Это достигнуто путём сочета- 
ния при построении модели АД двух подходов: схемотехнического для электрической части и ме- 
тода структурных схем [5] для магнитной и механической частей. 

В качестве датчиков и для ввода сигналов в составе модели АД использованы элементы 
типа ИНУН, ИНУТ, ИТУН. В этих сокращениях: И — источник, У — управляемый, Н — напряжение, 
Т — ток. Например, ИНУТ — источник напряжения, управляемый током. Графическое изображе- 
ние компьютерной модели фазы А обмоток статора и ротора АД, составленной по (1) — (5), пока- 


зано на рис. 1, где [.„ — индуктивность рассеяния фазы обмотки статора, [„ — индуктивность 


рассеяния фазы обмотки ротора. Традиционно параметры ротора приведены к обмотке статора. 
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ИНУТ! выполняет роль датчика тока фазы статора /.,, аналогично ИНУТ2 — роль датчика тока 
фазы статора /,. ИТУН1 служит для ввода в схему тока /, — приведённого тока фазы ротора. 
ИНУН1 вводит в цепь намагничивания фазы напряжение и,,. Аналогичную функцию выполняет 


ИНУНЗ — вводит в цепь обмотки ротора разность ЭДС ветви намагничивания и ЭДС вращения, то 
есть величину падения напряжения. ИНУН2 является дифференциальным датчиком разности по- 


тенциалов на индуктивности Г . Последующая обработка этого сигнала напряжения путём деле- 


а [е 
ния на значение /,„, даёт производную ет тока статора фазы. С помощью ИНУН4 выполняется 


[в] 
аналогичная операция с целью выделения сигнала рт | 


В компьютерной модели, изображённой на рис. 1, использованы сигналы ЭДС ветви на- 


магничивания (12) и ЭДС вращения (13) фаз ротора. 


@а = Мао 
ов = Иов, (12) 
е., =-И. 

лв 55 (4, -Ч, ), /\3; 

ев = (4, -4,), /\3; (13) 


е, = (4. -4,), /3. 


Активные сопротивления указанные на рис. 1, имеют значения порядка 10 МОм. Они 


Горел, 7 
введены в служебных целях. Практически не влияя на численные результаты расчёта, указанные 
сопротивления стабилизируют ход решения (моделирования) за счёт сохранения цепи тока (фи- 
зический смысл — путь протекания токов утечки). Это особенно важно при дискретном изменении 
сопротивления цепей питания АД, например при электрическом питании от полупроводникового 
преобразователя или при обрыве фазы. 

Адекватность динамики разработанной здесь компьютерной модели АД подтверждена хо- 
рошим совпадением результатов динамического моделирования на эталонных моделях АД (с ко- 
роткозамкнутым ротором) с номинальными мощностями 1,1 кВт и 4 кВт [1, 6]. Адекватность же 
статических режимов работы этой модели подтверждена многократными вычислительными экс- 
периментами для известных АД с короткозамкнутым ротором типов НВА-55 (номинальная мощ- 
ность 55 кВт) и АЖВ250М2ЗУХЛ2 (номинальная мощность 110 кВт) [7, 8]. Результаты этих вычис- 
лений удовлетворительно совпадают с данными расчётных записок и протоколов испытаний ОАО 
«ВЭлНИИ» (г. Новочеркасск) и ОАО «НИПТИЭМ» (г. Владимир). В настоящей работе описаны ре- 
зультаты компьютерного моделирования АД АЖВ250М23ЗУХЛ2 с параметрами и характеристиками, 


приведёнными в табл. 1. 


87 


Технические науки 












0 Ибо, > 
И НУТА и НУН 1 — ааа 






Рис. 1. Графическое изображение компьютерной модели фазы А обмоток статора и ротора АД 


Математическая и компьютерная модели ГМ. За основу модели ГМ взята физическая и ма- 
тематическая модель ГМ постоянного заполнения, описанная в [9]. Вращающий момент на валу 


привода насоса ГМ 
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2 
а гм > р (оса > 62 урбины "1 )о, (14) 
где О— расход рабочей жидкости в ГМ, м3/с, равный 
4 
Вы 90096,26 „А (=). 
2 Й 
Г. 15 
мя 1-6 и-9) (15) 
Г, 
= — относительное скольжение рабочих колёс ГМ, равное 
= = 1 р [бы ИО ); (16) 
ЛЕЁ (=) — опытная статическая характеристика ГМ, с”/м. 
Таблица 1 
Параметры и характеристики АД АЖВ250М23ЗУХЛ2, использованные при 
моделировании 
Обозначение Традиционное 
параметра обозначение 
№ Числовое Единица 
Наименование параметра в системе параметра в 
п/п значение измерения 
координат Т-образной схеме 
а, В, у замещения АД 
1 Количество пар полюсов р р 1 = 
Активное сопротивление фазы 
2 г П 0,0287 Ом 
обмотки статора 
Индуктивность рассеяния фазы 
3 [5 [а 0,000546 Гн 
обмотки статора 
4 Главная индуктивность фазы ГИ ГИ 0,021743 Гн 
Активное сопротивление цепи 
5 Г Г 137,051 Ом 
намагничивания фазы 
Активное сопротивление фазы 
6 обмотки ротора, приведённое к | г. [2 0,013 Ом 
статору 
Индуктивность рассеяния фазы 
7 обмотки ротора, приведённая к | [+ ь 0,000355 Гн 
статору 
8 Момент инерции ротора 3. 3. 0,484 кг"м? 
Номинальное линейное 
9 ет Ил лном 380 В 
напряжение статора 
Частота вращения в режиме иде- 
10 020 (2 ю 314,159 рад/с 
ального холостого хода 
Механические потери в режиме 
11 Риею Риею 740 Вт 
идеального холостого хода 
Описание других составляющих, участвующих в выражениях (14) — (18), дано ниже в 
табл. 2. 
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Уравнение движения для вала турбины ГМ имеет вид: 


[610 К 


турбины __ насоса ПМ `` сопр.Т 


(©) 


турбины = 


нагрузки 





СЕ У 


турбины 


+ 7кт 


+1 


нагрузки 


(17) 


В числителе правой части (17): первое слагаемое — вращающий момент на валу привода 


насоса ГМ (см. (14)); второе слагаемое — момент сил сопротивления вращению валов турбины ГМ 


и нагрузки; третье слагаемое — вращающий момент на силовом валу механизма внешней нагрузки. 


Таблица 2 


Значения параметров, использованные при моделировании силового 


электрогидродинамического привода 




































































№ Обозначение Числовое Единица 
Наименование параметра 
п/п параметра значение измерения 
1 Передаточное отношение редуктора р 1/1,2 - 
2 Момент инерции колеса насоса насоса 0,5 кг-м? 
3 Момент инерции жидкости насоса УЗкн 0,033 кг-м? 
Коэффициент, учитывающий увеличение момента 
4 инерции ротора АД за счёт приведённого момента | б 1,2 - 
инерции редуктора 
5 Момент инерции колеса турбины “Лурбины 0,5 кг-м? 
6 Момент инерции жидкости турбины Зет 0,033 кг-м? 
7 Момент инерции нагрузки нагрузки 2,4 кг м? 
Коэффициент  пропорциональности между частотой 
вращения вала АД в, и моментом сил сопротивления 
8 е Ксопр.ы 0,0095 = 
вращению валов насоса ГМ и редуктора, приведённым к 
валу ротора АД 
Коэффициент пропорциональности между частотой 
9 вращения вала турбины турбины и моментом сил сопро- | Копрт 0,0095 = 
тивления вращению валов турбины ГМ и нагрузки 
10 Ускорение свободного падения 9 9,81 м/с? 
11 Плотность рабочей жидкости ГМ р 850 кг/м? 
12 Активный диаметр ГМ р 0,363 м 
Малый радиус расчётного круга циркуляции (на средней 
13 й . п 0,054 м 
струйке) колёс ГМ 
Боьшой радиус расчётного круга циркуляции (на средней 
14 .. ы Б 0,154 м 
струйке) колёс ГМ 
Ширина кольца выходного сечения канала ведущего 
15 Ь 0,054 м 
колеса ГМ 
В (7) момент инерции / является суммой составляющих: 
,2 
о И (18) 


Для оценки (в первом приближении) временных характеристик быстропротекающих дина- 


мических процессов (например, при ударных нагрузках, вызывающих отрицательные угловые ус- 


90 





Вестник ДГТУ. 2012. № 1 (62), вып. 2 








корения порядка 2100 1/с?) последовательно с функцией Л =Ё (=), характеризующей статический 
коэффициент момента ГМ, установленный опытным путём, включено апериодическое звено 
К, (Е) ‚ учитывающее запаздывание во времени Ё реакции силовой гидродинамической передачи 
на скачок нагрузки, в виде: 
К, (+) =1-е 5. (19) 
Здесь &: — момент ударного нагружения ведомого вала ГМ, с; а — коэффициент, учиты- 
вающий динамику реакции ГМ на нагрузку, 1/с. Нами принято а? =2100 1/с? согласно рекоменда- 
циям [10]. 
М 


ном 





[4 





прямо 


Возведение модуля 
входного сигнала в 
степень 2 


© турбины 


> 


реверс 


Знак > 
входного сигнала 
нагрузки 


Рис. 2. Графическое изображение компьютерной модели способа формирования реактивного момента вязкого трения 


Чтобы придать разработанной здесь модели гидродинамической передачи свойства ГМ со 
статическим самоопоражниванием (тяговой) [11] с ограничением вращающего момента, переда- 
ваемого ГМ, удобно в компьютерной модели использовать нелинейность типа «насыщение». 

Для придания свойств ГМ с динамическим самоопоражниванием (предельной) [11], когда 
опоражнивание или заполнение рабочего объёма ГМ жидкостью начинается при переходе через 
критическое значение частоты вращения насоса ГМ ®насосакр, МОЖНО воспользоваться переходом с 


одной характеристики Л =Ё (=) на другую в зависимости от степени заполнения ГМ жидкостью. 
Принимая (для упрощения), что характеристики Л =#(=) при любом заполнении ГМ подобны по 
форме, можно домножать мгновенное значение функции Л =Ё (=), полученное на основной ха- 


рактеристике (при наибольшем коэффициенте заполнения ГМ жидкостью), на текущее значение 


коэффициента заполнения. Плавное изменение последнего при переходе через н„сосэкр ПОЛучим, 
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полагая процесс опоражнивания или заполнения жидкостью апериодическим с ПОСТОЯННОЙ вре- 


мени 0,2—0,3 с. Для этого воспользуемся апериодическим звеном типа ^, (Е) в виде: 


К (Е) =1-е*%). (20) 
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Рис. 3. Результаты компьютерного моделирования динамики привода для ГМ со статическим самоопоражниванием 
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Рис. 4. Динамическая характеристика Л =! (=) ‚ полученная в результате моделирования привода с ГМ со статическим 


самоопоражниванием при параметрах согласно табл. 1 и2 
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Здесь &› — момент достижения значения @несосакр Частоты вращения насоса ГМ, с; В — ко- 
эффициент, учитывающий динамику опоражнивания или заполнения рабочего объёма ГМ жидко- 


стью, 1/с. Нами принято 1/В =0,25 с. 


При моделировании /М/,„„рузки В случае ГМ со статическим самоопоражниванием использо- 

вался кусочно-постоянный сигнал, изменяющийся ступенчато (холостой ход, нагрузка, перегруз- 
ка). В случае ГМ с динамическим самоопоражниванием /„„ру»„ формировался как сумма реактив- 
ного момента вязкого трения (вентиляторного) и момента постоянной величины, действующего 
только в начальной стадии пуска, то есть до некоторой малой частоты вращения (момент трога- 
ния). На рис. 2 показан способ формирования реактивного момента вязкого трения, использован- 
ный в описываемой компьютерной модели. 
Результаты компьютерного моделирования. Результаты моделирования динамики силового 
привода показаны для ГМ со статическим самоопоражниванием на рис. 3. Здесь кривая 1 — час- 
тота вращения вала турбины @урвинь, Шс. Кривая 2 — частота вращения вала насоса ГМ ‘насоса, 
1/с. Последовательно показаны режимы разгона и работы на холостом ходу, работы под нагруз- 
кой, особенности работы в условиях перегрузки на валу турбины. 


А 
| 
Разгон АДи 
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Рис. 5. Результаты компьютерного моделирования динамики привода при ГМ с динамическим самоопоражниванием 


Для ГМ с динамическим самоопоражниванием результаты моделирования приведены на 
рис. 5, где график 1 — частота вращения турбины @хурвины, График 2 — частота вращения насоса 
ГМ @насос, График 3 — момент нагрузки на валу турбины М„„рузк, график 4 — момент на валу на- 
соса ГМ Мрнасосагм. Сценарий в этом случае следующий: пуск АД производится при незаполненном 
жидкостью рабочем объёме ГМ. Процесс заполнения ГМ рабочей жидкостью начинается при дос- 


тижении валом насоса частоты вращения 90 % от номинальной, то есть при достижении @насосакр 
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(это начало разгона нагрузки и турбины). Привод выходит на устойчивый установившийся режим 


работы. В момент времени Е, =10 с происходит отключение питания АД. Вследствие этого при- 


вод переходит в режим «выбега». При снижении частоты вращения вала насоса до 90 % от номи- 
нальной начинается процесс самоопоражнивания ГМ, по окончании которого ротор АД и нагрузка 
разъединяются, совершая далее «выбег» независимо друг от друга. 

Выводы. Результаты компьютерного моделирования электрогидродинамического привода адек- 
ватно отражают сущность реальных физических процессов [12], происходящих на динамических и 
статических режимах асинхронного электропривода с ГМ. Использование ГМ с динамическим са- 
моопоражниванием позволяет осуществить присоединение нагрузки к уже разогнавшемуся АД, 
существенно облегчая процесс его пуска [12], рис. 5. 
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